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5-H 1 ABSTRACT 

> 

Context. The space mission Kepler provides us with long and uninterrupted photometric time series of red giants. We 
are now able to probe the rotational behaviour in their deep interiors using the observations of mixed modes. 
Aims. We aim to measure the rotational splittings in red giants and to derive scaling relations for rotation related to 
seismic and fundamental stellar parameters. 

Methods. We have developed a dedicated method for automated measurements of the rotational splittings in a large 
number of red giants. Ensemble asteroseismology, namely the examination of a large number of red giants at different 
stages of their evolution, allows us to derive global information on stellar evolution. 
Results. We have measured rotational splittings in a sample of about 300 red giants. We have also shown that these 
splittings are dominated by the core rotation. Under the assumption that a linear analysis can provide the rotational 
■ splitting, we observe a small increase of the core rotation of stars ascending the red giant branch. Alternatively, an 

important slow down is observed for red-clump stars compared to the red giant branch. We also show that, at fixed 
| stellar radius, the specific angular momentum increases with increasing stellar mass. 

Conclusions. Ensemble asteroseismology indicates what has been indirectly suspected for a while: our interpretation of 
the observed rotational splittings leads to the conclusion that the mean core rotation significantly slows down during 
the red giant phase. The slow-down occurs in the last stages of the red giant branch. This spinning down explains, for 
instance, the long rotation periods measured in white dwarfs. 

<-h \ Key words. Stars: oscillations - Stars: interiors - Stars: rotation - Stars: late-type 

1. Introduction ments of the surfa ce rotation are giv en by the mea- 
sure of v sini (e.g. ICarnev et all 120081 ). The slow rota- 

The internal structure of red giants bears the history of tion rate in i ow _ mass wn ite dwarfs (e.g. iKawaler et all 15 

their evolution. They are therefore seen as key for the un- rjgggj) SU g ge sts a spinning down of the rotation during the 

derstanding of stellar evolution. They are expected to have red giant branch ( RGB ) phase In addition, 3-D simula- 

5 a rapidly rotating core and a slowly rotating envelope (e.g. tions show n 0n - r i g id rotation in th e convective envelope of 

I Sills fc Pinsonneault! | 2000|) , as a result of internal angular red giants flBrun fc Palaciosl l2Pjp^ . Different mechanisms 

momentum distribution. Indirect indications of the internal fo r spinning down the co re have been investigated (e.g. 20 

angular momentum are given by surface-abundance anoma- ICharbonnel fc Talonl l2005h . Rotationally-induced mixing, 

lies resulting from the action of internal transport processes amid other angu i ar momentum transport mechanisms, is 

10 and from the redistribution of angular momentum and st ill poorly understood but is known to take place in stellar 



chemical elements (|Zahn 1992; Ta on fc Charbonnel||2008j; interiors. Therefore, a direct measurement of rotation in- 

iMaedaJ l2009t | Canto Martins et al. | 12011D . Direct measure- side re d giants would give us an unprecedented opportunity 25 

to perform a leap forward on our understanding of angular 
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Fig. 1. Echelle diagram of a typical RGB star (KIC 
6144777) as a function of v/ Av — (n p + e). The radial order 
n p is indicated on the y-axis. Radial modes (highlighted in 
red) are centered on 0, quadrupole modes (highlighted in 
green), near —0.12 (with a radial order n p — 1), and I — 3 
modes, sometimes observed, (highlighted in light blue) near 
0.20. Dipole mixed modes are identified with the frequency 
given by the asymptotic relation of mixed modes, in /iHz. 
The fit is based on peaks showing a height larger than eight 
times the mean background value (grey dashed lines). 



momentum transport in stellar interiors (e.g. lLaearde et all 
l2012HEggenberger et aLll2012ft . 

This is becoming possible with seismology, which pro- 
vides us with direct access to measure the internal rotation 
profile, as shown bv lBeck et al.l (|2012t ) and lDeheuvels et all 
(|2012ah . They demonstrate that it is possible to measure 
the core rotation thanks to the precision of asteroseismic 
results derived from the photometric light c urves of red gi- 
ants provided by the NASA Kepler mission (jBorucki et al.l 
I2010T ). Previous work on field star asteroseismic surveys 
shows that is possible to define evoluti onary sequences (e.g. 
iMielio et al1l2009t lluber et al.ll201lft . So, with ensemble 
asteroseismology, we aim to get measurements in a large 
enough sample of red giants to explore the variation of the 
internal rotation with evolution. 

After a decade of relatively uncertain ground-based 
measurements, CoRoT has unambiguous ly revealed that 
red giants show solar-like oscillations (|De Ridder et al.l 
l2009t). Important sc aling relations have then been shown 
(Hekkcr et al.l l2009f) for deriving crucial information on 
the s tellar mass and radiu s from global seismic parame- 
ters (Kallinge r et alJ 12010). Specific features of solar-like 
oscillations in red giants have been characterized (e.g . 
IBedding et aT1l2010l: iMosser et aLkOlOt iHuber et al.ll201oh . 
Mixed modes, which correspond to the coupling of grav- 
ity waves in the r adiative core region and pressure waves 
in th e envelope (jDziembowski et all 120011 : iDupret et al.1 
l2009f ). have been detected in red giants. They were first 



reported by IBedding et al.l (l2010t). Their period spacings 55 
were measured bv lBeck et alJ (|2011l) . IBedding etafl (|2011l) 
and Mosser et al. (2011a) have shown the capability of 
these modes to measure the evolutionary status of red 
giants. Mixed modes can be divided into two categories, 
namely gravity-dominated mixed modes (hereafter called 60 
g-m modes) which have large amplitudes in the core, and, 
in contrast, pressure-dominated mixed modes (p-m modes). 
The frequencies of the p-m modes are very close to the the- 
oretical pure p mode frequencies; they however appear to 
have a significant g component, and are therefore sensitive 65 
also to the core conditions. We analyse in this work mostly 
stars showing a rich mixed-mode spectrum. 

Observations and data are presented in Section [2] The 
observed properties of the rotational splittings are de- 
scribed in Section[31 In Section^ we derive scaling relations 70 
governing the rotational splitting, independent of any mod- 
eling, but based on the observational evidence of a much 
higher rotation rate in the stellar core. The way the core 
rotation is related to the measured rotational splitting is 
quantified in Section [SJ In Section [51 we then derive unique 75 
information on the core rotation in red giants and probe 
their internal angular momentum and its evolution. 



2. Data 

2.1. 25-month long observation 

The red giant stars analyzed in this wo rk have already 80 
been presented (e.g. iHekker et al.l I2011L and references 
therein). We now benefit from longer time series. All red 
giants observed up to Kepler's quarter Q10 have been an- 
alyzed. Original ligh t curves were processed according to 
Uenkins et al. (I2010D and cor rected according to the pro- 85 
cedure of Garcia et al.1 (|201ll) . The Fourier analysis of the 
868-day long time series provides a frequency resolution of 
about 11.5 nHz. Due to the characteristics of the measure- 
ment, we have in principle access to rotation periods up 
to the observation duration. For a few RGB stars showing 90 
large rotational splittings, a reprocessed set of Q0 to Qll 
Kepler data has been used, based on the extraction of the 
stellar fluxes using new custom masks from the recently re- 
leased pixel-data information (Bloemen et al. 2012). These 
new light curves w ere then corrected applying the algo- 95 
rithms developed bv lGarcia et al.l(|201lD but using a refined 
new automatic procedure (Mathur et al., in preparation). 
In practice, we measure rotational splittings with periods 
in the range 8 - 280 days. 



2.2. Mixed-mode pattern 

The complete identification of the red giant pressure oscil- 
lation pattern is given by the description o f the so -called 
universal oscillation pattern (|Mosser "etaLl l20lTbh . This 
method alleviates any problem of mode identification. The 
whole frequency pattern of pure p modes is approximated 
by: 
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where Av is the mean large separation measured in a 
broad frequency range around the frequency f max of max- 
imum power, rip is the p-mode radial order, £ is the an- 110 
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gular degree, e is the phase offset, doe accounts for the so- 
called small separation, a is a small constant, and n max = 
v max /Au. The parameters e, doe and a are considered as 
a function of the large separation; an updated fit of e i s 
used, comparable to the expression of lCorsaro et~aT1(l20Tl . 
The parameter a represe nts the second -order term of the 
asymptotic development (jTassoull [19 801) and accounts for 
the mean curvature of the radial mode oscillation pattern 
(IMosser et al.l|2012a|K It was considered as a constant by 
M osser et al.l (|2011bh . Here, with much longer time series 
and large separations observed up to 20 /iHz, we prefer to 
use the fit 



= 0.015 At/ 



-0.32 



(2) 



with Av in y«Hz. This relation is derived from the detailed 
analysis of the radial modes with the method presented by 
iMosserl (|2010H . 

The recent analysis of iKallinger et all ()2012l ) has shown 
that this development is valid, independent of the stellar 
evolutionary status, under the condition that the determi- 
nation of the large separation is global and not local. Their 
Figure 6 illustrates the important curvature of the ridges 
depicte d by the quadratic te rm of Eq. (JlJ . The comparative 
work by Hekkcr et al. (201 2J) sho ws that the extra hypothe- 
sis used bv lMosser et al.l (|2011bl ) , expressed by the function 
e(Ai/), helps to obtain very precise results. Using this con- 
straint requires fitting each oscillation spectrum in a large 
frequency range around ^ ma x and taking into account the 
mean curvature of the p-mode oscillation pattern. 

Equation §T§ holds precisely for radial modes and gives a 
proxy for the pure pressure dipole (£ = 1) modes. However, 
due to the significant coupling of pressure waves with grav- 
ity waves in the inner radiative region, the dipole oscillation 
pattern is dominated by mixed modes located around the 
position of the pure p mode. The frequency interval between 
the individual mixed modes is not constant and is deter- 
mined according to the method presented in iMosser et al.l 
(2011a). For stars showing a large number of g-m modes, 
the oscillation pattern can be p recisely descr i bed b y an 
asymptotic relation presented bv lGoupil et al.l (2012), fol- 
low ing the ideas origin ally developed by Shibahashil (1979) 
and Unno et al. (Il989h. We apply this method, as explained 
m M osser et al.l (l2012dV in order to locate precisely the 
dipole mixed modes. 

The echelle diagram of Fig.Q]shows the identification of 
the radial modes provided by Eq. ([1]), with the mode cur- 
vature, and the location of the mixed t = 1, m = modes 
defined by the asymptotic relation. The remaining shifts 
between the actual and expected peak s positions, due fo r 
instance to a sharp structure variation (Mig lio et alJl201Cft , 
are small compared to the mixed mode spacings. Hence, 
they do not hamper the mode identification. 



3. Measuring the rotational splittings 



From the analysis of dipole mixed modes. iBeck et aD (|2012l ) 
showed differential rotation in red giants. In this Section, 
we propose a model for describing the rotational splittings 
of dipole mixed modes and an automated method for mea- 
suring them. 
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Fig. 2. Empirical rotation profile TZ n for dipole mixed 
modes associated to the pure p mode of radial order n p , 
as a function of the reduced frequency v / Ai> — (n p + e). 



3.1. Modulation of rotational splittings 

Because of differential rotation, a simple interpretation of 
rotational splittings in terms of mean rotation of the star 
is inadequate. The situation is made even more compli- 
cated by the fact that rotational splittings are measured 
for dipole mixed modes. Indeed, the mixed nature of the 
m odes varies as a fun ction of mode frequency, as shown 
by iDupret et al.l (120091) . As a result of those two interlaced 
effects. IMosser et al.l (|2012d ) have shown that rotational 
splittings are modulated in frequency, with a period of Av. 
The rotational splitting can be written 
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'split 



''n,l,m - v n,l = m H{v) 5v rot , 



(3) 



as: 
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•Vl,™ = v n,l + m ft(fn,l,m) ^r, 
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175 



where m is the azimuthal order and Sis TOt is the maximum 
value observed for g-m modes. Locally, aroun d the position 
of a pu re dipole p mode of radial order n p , IMosser et al.l 
(2012cj) have shown that 1Z can be empirically expressed 



180 



(4) 185 



for a mixed mode with frequency v associated with the pres- 
sure radial order n p . The development, solely derived from 
the observed splittings, has currently no theoretical basis: 
the Lorentzian form has been chosen similar to the ob- 
served variation of the mixed-mode sp acing with frequency 
(IBeck et al.H2011l:lBedding et al.ll201lh . 

The form of TZ (Fig. [5]) implies that the splittings are not 
symme tric, consistent with the findings of iDeheuvels et ail 
(|2012bl) . For multiplets near p-m modes, the closest compo- 
nent to the theoretical pure p mode has the smallest split- 
ting. The asymmetry is expressed when considering Eq. 
as an implicit equation, the splitting of the frequency v n ^ m 
compared to the non-rotating reference v n ^\ depending on 



190 



195 



(5) 200 



It assumes, as observed but with the limitation of the ob- 
served frequency resolution, that the modulation of the 
splitting is the same for all radial orders. In fact, the fits of 
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the splittings, obtained by supposing that the frequency of 
205 the m = component of the dipole mixed modes is given 
by the asymptotic development, show that the terms A and 
/3 are independent of the frequency. Furthermore, we have 
verified that the values of A and j3 are always very close to 
0.5 and 0.08, respectively, and hence do not depend on the 
210 evolutionary stage. 

Due to the large volume of data, we need an auto- 
mated method for deriving the rotational splittings, sim- 
ilar to the determi nation of the globa l seismic parameter s 
Av and ^ max (e.g. iHekker et aljfeoill IVerner et al.ll201ll) . 
215 This method has to cope with an interweaving of rotational 
splittings and mixed-mode spacings. 



3.2. EACF analysis 

We can use the envelope autocorrelatio n function (EACF) 
method (Mo sser fc Appourchauxl [20091 ) to measure the ro- 

220 tational splittings of non-radial modes. This method was 
developed to measure the frequency spacing correspond- 
ing to the large separation of solar-like oscillations. With 
narrow filters centered on the expected 1 = 1 pressure 
mode s, it gives the spaci ng due to the mixed-mode pat- 

225 tern osser et all 1201 lal ). With ultra -narrow filters cen- 
tered on each individual mixed mode, it proves to be able 
to measure the rotational splittings. The central positions 
of the filters are chosen within the frequency ranges where 
mixed modes are expected. Each ra nge is wider than ha lf 

230 the large separation Av, as shown bv lMosser et all (|2012bh . 
Different central positions of the filters are tested within 
this range, independent of the mixed-mode positions. The 
widths of the filters have to be narrow enough to select only 
one mixed mode, in order to avoid confusion between the 

235 g-mode spacing and the rotational splitting (Fig. [3]). They 
are varied, in order to test a large range of splittings. We 
have tested five different widths, of the order of 1 //Hz or 
less, which encompass the observed values. 

In the EACF method, the signature of any comb-like 

240 structure in the spectrum, here the rotational splitting, is 
provided by the highest peak in the spectrum of the win- 
dowed spectrum. The signal is normalized in such a way 
that the mean white-noise level is 1. The reliability of the 
detection is derived from an HO test. Since we use ultra- 

245 narrow filters, only high signal-to-noise time series can be 
analyzed. 
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v GuHz) 

Fig. 3. Zoom on the oscillation spectrum of the target KIC 
10777816. Different narrow filters centered in the 1 = 1 
mixed mode range, indicated with different line styles, al- 
low us to measure a local rotational splitting in each filter. 
For clarity, only those filters centered on possible multiplets 
have been represented. 
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Fig. 4. Comparison of the splittings measured with the 
EACF automated method to the splittings measured with 
the fit of the mixed modes (Eq. ([I])). Red asterisks cor- 
respond to rotational splittings too large to be accurately 
measured in an automated way and hence excluded from 
the sample. 



For each filter width, multiplets were searched in seven 
radial orders around f max , at ten different positions per ra- 
dial order. As a consequence, seventy possible signatures 
of the rotational splittings are analysed for a given filter. 250 
The criterion for a positive detection is the measurement 
of similar splittings for at least three positions centered 
on different dipole mixed modes corresponding to differ- 
ent radial orders n D , with a signature of the EAC F above 
the threshold level (jMosser fc Appourchauxl [20091 ) . In or- 255 
der to account for the modulation of TZ towards p-m modes 
(Fig. [5]), we allow relative differences of 1/3 between the 
individual measurements. This threshold level has been de- 
termined empirically. It takes into account the fact that 
splittings are measured mainly in the wings of the function 260 
TZ: the splittings of p-m modes are not easily measured, due 
to their shorter lifetime, whereas the amplitudes of the g-m 
modes located far from the p-m modes are too low to give 
a reliable signature. 

3.3. Performance 265 

All results found by the automated method have been ver- 
ified individually by visually comparing the spectrum with 
itself after a shift corresponding to the rotational splitting. 
The method appeared to be dominated by g-m modes since 
they have narrower widths and are more numerous than 270 
p-m modes. Hence, they contribute much more efficiently 
to the EACF signature, as shown by the careful examina- 
tion of numerous individual spectra. Therefore, the method 
mainly gives the rotational signature of the core. This is also 
clear from the close examination of the extracted splittings. 275 

For stars with a low signal-to-noise ratio oscillation 
spectrum, generally faint stars or stars with low i/ max , we 
sometimes obtain spurious results. These results are clearly 
caused by the stochastic nature of the oscillation excitation, 
which occasionally resembles the complex pattern of mul- 280 
tiple peaks shown by the dipole mode, and were discarded. 
Finally, depending on the stellar inclination i, the multi- 
plets have two components (when sini is close to 1), three 
components (intermediate i values) , or only one component 
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Fig. 5. Zoom on the rotational splittings of the mixed 
modes corresponding to the radial orders n p — 8 — > 11 
in the giant KIC 9574650, in an echelle diagram as a func- 
tion of the reduced frequency v/ Av — (n p + s). All triplets 
(m = —1,0, 1) are identified with the function 1Z defined 
by Eq. ([4]). At low frequency, multiplets are overlapping. 



285 (low i). In this latter case, measuring the splitting is not 
possible. The correct identification of the multiplets has 
been possible in most of the cases, allowing us to remove 
the confusion between 5z^ sp iit and 2 <5^ sp iit- 

Possible confusion between rotational splittings and 

290 mixed-mode spacings has been investigated. Such a con- 
fusion is usually eliminated by limiting the width of the 
filter. However, RGB stars may show a rotational splitting 
very close to the g-mode spacing. We therefore explored 
the frequency domain where the solutions are ambiguous. 

295 Ambiguous detections are identified by the fact that, even 
if the rotational splitting and the mixed-mode spacing are 
both modulated in frequency, with a period Av, their sig- 
natures are different. On the one hand, the g-mode fre- 
quency spacing varies as v 2 since it approximately corre- 

300 sponds to a regular spacing in period (jBedding et al.ll201lt 
iMosser et al.ll2012d ); on the other hand, the modulation 1Z 
is the same all along the spectrum (Eq. 

As expected, the automated method fails when the ro- 
tational splittings are larger than half the mixed-mode 

305 spacings. This occurs for giants ascending the RGB, with 
Av < 12 //Hz. For these stars, a dedicated method, pre- 
sented in the next paragraph, is needed to disentangle the 
rotational splittings from the mixed-mode spacings. We fi- 
nally measured reliable rotational splittings in 265 red gi- 

310 ants with the automated method, in the clump and in the 
early stages of the RGB. The combined effect of the tiny 
width of the filter and of the limited frequency resolution 
makes the method much more precise at high frequency 
than at low frequency. The relative precision is of about 

315 5% in 5i/ ro t for stars with Av — 15/iHz and 25% when 
Av = 5 /iHz. 



ation in red giants 

Table 1. Positive detections as a function of Av 



Av range Reference Detections 



( M Hz) 


(a) 


(b) 


(c) 


(d) 


(e) 


3 


323 


323 


49 





0.0% 


3 4 


257 


250 


76 


37 


48.7% 


4 5 


412 


387 


167 


146 


87.4 % 


5 6 


97 


86 


82 


21 


25.6 % 


6 8 


92 


81 


75 


31 


41.3% 


8 12 


68 


51 


49 


48 


98.0 % 


12 20 


50 


38 


36 


30 


83.3 % 


all 


1299 


1216 


534 


313 





(a) : total number of red giants in a given frequency range. 

(b) : same as case (o), but spectra with depressed dipole mixed 
modes are excluded. 

(c) : same as case (6), but spectra with an EACF signature less 
than 100 are also excluded. 

(d) : number of positive detections in each frequency range. 

(e) : percentage of positive detections with respect to reference 
(c). 



3.4. Validation with direct measurements 

The values of the rotational splitting provided by the EACF 
method can be compared to a direct fit of the modulation of 
the splitting based on Eqs. ((U) and ([3]). A tutorial for fitting 320 
the rotational multiplets is given in Appendix [X] different 
fits are shown, in order to illustrate cases with rotational 
splittings lower, comparable or larger than the mixed-mode 
spacings. Using such fits, we measured the rotational split- 
ting Sv Iot in 102 red giants. Among them, 54 had already a 325 
rotational splitting measured with the EACF method. The 
remaining 48 stars had rotational splittings too large com- 
pared to the mixed modes spacings to be measurable with 
the EACF method. As a result, the total number of red 
giant stars with splittings measured with one or the other 330 
method is 313 (=265+102-54). The analysis of the detec- 
tion and non-detection as a function of the large separation 
is summarized in Table [TJ Non-detection at large Av oc- 
cu rs for RGB stars wit h depressed mixed modes (identified 
bv lMosser et al.ll2012b[ ) or with very low inclination. In this 335 
latter case, the components m = ±1 are too low to allow the 
identification of multiplets. At small Av, the non-detection 
is explained by the limited frequency resolution, but also 
by the poorer qua lity of the spectra, as express ed by the 
EACF coefficient (|Mosser fc Appourchauxl 12009') . For Av 340 
in the ranges [4, 5/iHz] and [5, 6/iHz], the success of the 
detection depends mainly on the evolutionary status: small 
rotational splittings in clump stars with large mixed-mode 
spacings can be detected more easily than larger rotational 
splittings of RGB stars embedded in narrow spacings. We 345 
expect the number of positive detections to increase with 
prolonged observations: the frequency resolution will give 
access to rotational splittings in the upper RGB, and the 
better signal-to-noise ratio will allow the identification of 
tiny m = ±1 components at low stellar inclination. 350 

We have noted that the automated EACF approach pro- 
vides values that are 10 % smaller than Sv mt derived from 
the individual fits (Fig. [4]). This small correction is consis- 
tent with the different principles of the methods: the in- 
ferred value of Sv Iot from Eqs. ((3]) and (jH) is larger than 355 
any observed splitting, whereas the EACF method, even if 
dominated by g-m modes, also includes narrower multiplets 
in the vicinity of the p-m modes. This indicates that the 
10 % difference can be considered as a bias of the automated 
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360 method. For homogeneity, all splittings obtained with the 
EACF method have been multiplied by a factor 1.10. 

We have tested that the differential-rotation term 1Z of 
Eq. (j4]) holds for red giants in all stages of their evolu- 
tion. Residuals of the fit of the splittings are much smaller 

365 than 6v mt . We observe a coefficient A (Eq. (U])) of 0.5±0.1 
for an early RGB star such as KIC 7341231 studied by 
iDeheuvels et al.l (|2012al) . with Av = 28.9 /iHz, as well as 
for clump stars with Av ~ 4.0 /iHz. The parameter /3 is 
about 0.08 ± 0.015 for all giants, except a very small num- 

370 ber of exceptions, at the bottom of the RGB. Measuring a 
modulation profile 1Z almost independent of the evolution 
may indicate that its origin obeys to generic properties, 
similar for all red giants ascending the red giant branch 
(RGB). 

375 For clump stars and early RGB stars, we note that 
the rotational splittings are significantly smaller than the 
mixed-mode spacings. In these cases, the individual fits con- 
firm the automated measurement. However, when the rota- 
tional splittings are larger than the mixed-mode spacings, 

380 the EACF method was ineffective for extracting the cor- 
rect large splitting: it derived wrong values that result from 
a combination of the mixed-mode spacing and rotational 
splitting. Application of Eq. allowed us to provide the 
correct splittings even with large values comparable to the 

385 mixed-mode spacing. We stress that, in practice, it seems 
impossible to disentangle such multiplets from the mixed- 
mode pattern without the use of the asymptotic relation 
of mixed modes. We were able to investigate complicated 
cases, with rotational splittings significantly larger than 

390 the mixed-mode spacings, for RGB stars with Av down 
to 7 /iHz. Finally, even in the case where multiplets over- 
lap, that is zAi m +i,i.-i < ^n m ,i.+i, with n m the mixed-mode 
order, we do not observe any modification of the profile TZ 
(Fig. [5]). This indicates the absence of avoided crossings 

395 between dipole mixed modes with consecutive mixed-mode 
orders n m and different azimuthal orders to. In other words, 
this shows the absence of coupling between such modes. 

4. Scaling relations 

We have derived estimates of the stellar mass and ra- 
400 dius from the seismic global parameters and from the ef- 
fective temperatur es given by the Kepler Input Catalog 
(jBrown et al.l 1201 ll ). so that we can present the rotational 
splitting as a function of the stellar radius R (Fig. [7]). 
The evolutionary st atus of the stars was determined by 
405 iMosser et al.l (|2012ch . Among the helium-burning stars, we 
consider those with a mass greater than 1.8 M Q to belong 
to the secondary clump. The uncertainties arising from the 
effective temperatures and the imprecision of the scaling 
relations do not impact the following analysis. 
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Fig. 6. Rotational splitting 5 function of the large 

separation Av, in log-log scale. The dotted line indicates 
the frequency resolution. The dashed and dot-dashed lines 
represent the confusion limit with mixed m odes in RGB 
and clu mp stars, respectively, derived from IMosser et al.l 
(|2012ch . Crosses correspond to RGB stars, triangles to 
clump stars, and squares to secondary clump stars. The 
color code gives the mass estimated from the asteroseis- 
mic global parameters. The vertical bars indicate the mean 
error bars, as a function of the rotational splitting. 




Fig. 7. Rotation splitting S function of the aster- 

oseismic stellar radius, in log-log scale. Same symbols and 
color code as in Fig. [51 The dotted line indicates a split- 
ting varying as R~ 2 . The dashed (dot-dashed, triple-dot- 
dashed) line corresponds to the fit of RGB (clump, sec- 
ondary clump) splittings. 



410 4.1. Rotational splittings Sv YOt 

Rotational splittings Sv ro t are shown as a function of the 
large separation Av (Fig. E]). We see that the detection is 
difficult at low Av, due to the limited frequency resolution. 
This precludes an analysis of the core rotation in the high- 
415 luminosity RGB stars, but allows us to measure rotation in 
the clump stars. 

We first consider the RGB stars, indicated by crosses 
in Fig. El We note that the rotational splitting slightly de- 
creases when Av decreases, that is, when the star evolves 



on the RGB. For clump stars and secondary clump stars, 420 
splittings are much smaller. At this stage, ensemble aster- 
oseismology indicates either that the core rotation spins 
down, or that the splitting of g-m modes is not dominated 
by the core rotation. This last result is unlikely since g- 
m mode splittings are significantly larger than p-m mode 425 
split tings (Fig p> , co nsistent with the observations reported 
bv lBeck et all (|2012ft and IDeheuvels et al.l (j2012al) for four 
early RGB stars. 
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4.2. Scaling relations 

430 We have examined how the rotational splittings evolve with 
the stellar radius along the RGB (Fig. EJ). Only RGB stars 
with a mass in the range [1.2, 1.5 M©] were considered, to 
avoid a bias from the fact that high-mass stars are under- 
represented in the early stages of the RGB, whereas low- 

435 mass giants are under-represented in the later stages. With 
49 RGB stars in this case, we find 

Sv rot oc i?-0-5±o.3 ( RGB) _ ( 6) 

In the first stages of the RGB, the splittings <5^ rot show a 
slow decrease. Assuming the local conservation of angular 
440 momentum, such a decrease seems in contradiction with 
the core contraction: this has to be investigated. 

The same exercise can be done for the clump stars. 
The fit, conducted over a much broader range of mass 
(jMosser et al.ll2012d) . gives 

445 8v mt cx i?- 13±a4 (clump) or cx R- 1A±0A (2nd clump). (7) 

This indicates a different behaviour compared to RGB 
stars. We first note that the slopes are independent of the 
stellar mass. Closer to —2 than for RGB stars, they cer- 
tainly relate the influence of the stellar expansion. We also 
450 note that secondary clump stars, which are more massive, 
show larger splittings than clump stars. 

4.3. Intermediate conclusions 

From the analysis of the rotational splittings with the 
stellar radius, we note the weak decrease of RGB stars. 

455 According to the exponent of the fit reported by Eq. ([5]), 
this g-m mode splitting cannot be related to the sur- 
face rotation if it evolves at constant local angular mo- 
mentum. The large change in the rotation evolution from 
the RGB to the clump can be related to th e expansion 

460 of the non-degenerate he lium burning core (jlbenl fl97TL 
ISills fc Pinsonneaultl[2000l) . This increase of the core ra- 
dius is however limited and cannot explain the entire ob- 
served decrease of the rotational splittings, so that we are 
left with the most plausible conclusion that the strong de- 

465 crease of the rotational splitting is the signature of a signifi- 
cant transfer of internal angular momentum from the inner 
to the outer layers. This transfer should preferably occur at 
the tip of the RGB, out of reach with current Kepler ob- 
servations due to a limited frequency resolution. One could 

470 also imagine that the rotational splittings are sensitive to 
different layers, depending on the evolutionary status. This 
is investigated in the next Section, where we aim to inter- 
pret the signification of the observed splittings Sv ro t . 

5. Linking the rotational splittings to the core 
475 rotation 

To go a step further, we intend to qualitatively link the 
observed rotational splittings to the rotation inside the red 
giants. 



5.1. Linear rotational splittings and average rotation 

480 We assume that the rotation is slow enough that a 
first-order perturbation theory is sufficient to com- 
pute the rotational splittings. This yields t he fol- 
lowing expression for rotational splittings (|Ledouxl 



1951j Christensen -Dalsaaard & Frandsen 

Christensen-Dalsgaard fc Berthomieul 119911 : iGoupill 
iGoupil et al.ll201 

a f 1 r ( \ a 

ov mt , n ,e = / K-n,l\%) — — dx, 



1983 



2009; 485 



2tt 



(8) 



where x = r/R is the normalized radius and is the an- 
gular rotation (in rad/s). The rotational kernel JC n .i of the 
mode of radial order n and angular degree I takes the form 490 



fc n ,e — 



1 

I'll.! 



4 2 + (A - 1)£ _ 



px 



where I n / is the mode inertia 



Ird — 



& + A £ px 2 dx. 



(9) 



(10) 



The quantities entering these equations are the fluid verti- 
cal and horizontal displacement eigenfunctions, £ r , £h re- 
spectively; p is the density, and A = £(£ +1). 

The mode inertia is much larger for g-m modes that 
have a large amplitude in the inner cavity than for p-m 
modes. For a red giant, several mixed modes exist in the 
frequency vicinity of each radial mode. As a result, the vari- 
ation with frequency of the dipole mode inertia shows a reg- 
ular variation wit h a periodicity roughly eq u al to the large 



separ a tion (see LD zicmbowski et al.l 120011: iDupret et al.l 



2009; IChristensen-Dalsgaardl |20 1 It IGoupil et all 12012ft. 



This has been observationally confirmed (jMosser et al.l 
|2012(J| . The linear rotational splittings (Eq. ©) are found 
to follow closely the same behavior as the mode inertia for 
the same reason. However a significant amplification of the 
variation of the splittings with frequency compared to that 
of mode inertia sh ould exist when th e rotation is large in 
the central region (jOoupil et al.H2012h . 

Because the red giants are characterized by an inner 
dense region and an outer envelope, it is convenient to con- 
sider the rotational splittings as the sum of two contribu- 
tions 



Z7T 



<n K )c 



(ii) 



(^K)coro,«,£ = / O(z) K n jXx) dx, 

JO 

(fijOenv.rM = / ^( x ) ^n,t( X ) Ax - 



(12) 
(13) 



The core boundary x corc must be understood here as the 
limit where Q(x) K n ^(x) no longer contributes to the in- 
tegrant in (fiK)core, n,£ ■ Numerical c alculations show tha t 
a^coro remains the same for all modes (|Maraues et al.ll2012j ) . 
Equation (jlip is then equivalent to 



8lA 



rot ,n,£ 



2tt 



1 + a rot ) 



where a rot can be written 



«rot 



A" 



A' 



(14) 



(15) 



495 



500 



505 



510 



515 



where (OK) C ore,n,£ is the angular rotation, weighted by the 
kernel, averaged over the central layers enclosed within a 
radius r corc and {^lK)cnv,n,£ the angular rotation averaged 
over the layers above x core = r corc / R, 520 



525 
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535 



540 



545 



550 



555 



560 



with the definitions 



Kc 



K, 



env, 



Kn,t(x) dx, 



]C n> t(x) dx, 



530 and 



J Xcore fi(x) K n ,e(x) dx _ (f2 K ) c 



K,, 



Q(x) Kn t e(x) dx 



(16) 



(17) 



(18) 



(19) 



We then consider dipole g-m modes that have the largest 
rotational splittings ; they coincide with the largest inertia 
(Dup ret et al.1 12003 ). These modes are furthest away from 
the nominal pure pressure dipole modes and close to the 
radial modes. The splittings associated with these modes 
do not vary with frequency, i.e. from one radial mode to the 
other. Hence, we retrieve the observed splittings (Eqs. (j3|) 
and (HD): 



Si/ rot = max(di/ ro t,rc,£=iJ = max{dv sp i it 



(Si 



(20) 



As a consequence, we can drop the subscripts n, £ from now 
on. One expects, for red giants, a more rapid rotation rate 
in the inner layers hence (O) e nv/(^)core < 1 and <C 1 for 
very fast rotating cores (jGoupil et alJ 12012). For the g-m 
modes, numerical calculations show that K env / K cole -C 1 
(Fig. ED, hence a rot <C 1. Thus, from Eqs. JI3D and (fI5|): 



5v r . 



(n K ) 



<n>c 



2tt 



2tt 



(21) 



For these modes, the displacement is essentially horizontal 
in the core, therefore £ r <C £f We also have I ~ icorej so 
that, from Eqs. ([§]) and (|10l) . one can derive that the core 
kernel reduces to abou t 1/2, in agree ment with the Ledoux 
coefficient of g modes (|Ledouxlll95ltl . Finally, we have 



2 TV 



(22) 



Hence, rotational splittings of g-dominated modes provide 
a measure of the rotation averaged over the central region. 
The averaged rotation roughly corresponds to the value of 
the rotation at the radius where the mode amplitude of 
the horizontal displacement £h is maximum. This happens 
away from the center, in a core region where the rotation 
can have significantly decreased compared to the central 
rotation. In that case, the average rotation value gives a 
lower limit of the rotation of the very deep layers. If the 
rotation happens to be nearly solid in the central region, 
then the average rotation gives the rotation of the nearly 
uniformly rotating core. 

Keeping in mind these limits, <5^ ro t can be considered as 
a proxy of the mean rotation period of the core (T rot ) c , i.e. 



565 (T rot ) c = 



2tt 



1 



(fi) core ^ O^rot 

for dipole mixed modes. 



(23) 
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Fig. 8. Top: mode inertia for radial and dipole modes 
of a 1 Mq RGB star at the bump. T he evolution model 
has been calculated with CESAM (jMorell Il997t) and 
the oscillation frequencies w ere obtained with ADIPLS 
(|Christensen-Dalsgaardl |2011[ ). Bottom: normalized inte- 
grated rotational kernels of three dipole modes: M\ and 
M2 are g-m modes, whereas M3 is a p-m mode. The ver- 
tical triple-dot-dashed line indicates the mean location of 
the hydrogen-burning shell, and the vertical dashed line in- 
dicates the base of the convective envelope. 



5.2. Information from the kernels 



With 
1991 



the 
and 



CESAM code for 
the ADIPLS code 



stellar evolution (jMorell 
for adiabatic oscillations 
(Christcnscn-Dalsgaard] 1201 ll) . we estimate the kernels 570 
in red giant models at different evolutionary states. 
Figure [S] shows the normalized integrated rotational kernels 
J Q r /C(r) dr/ Jq ? X(r) dr derived for a p-m mode and two g-m 
modes in an RGB star at the bump, with At/ ~ 5 //Hz. We 
verify that the kernels in g-m modes are dominated by the 575 
core, since the normalized integrated kernels reach a value 
larger than 0.95 at the core boundary. 

These integrated kernels allow us to derive a refined es- 
timate of the core rotation. With a 2-layer model presented 
in Appendix IB. 11 a small correction to Eq. (|2^|) can be in- 580 
troduced by a factor rj (Eq. (|B.3p ) 



(T TO t)c 



(24) 
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585 



590 



595 



The correction factor r/, slightly larger than unity, accounts 
for the propagation of the oscillation in the envelope. This 
correction is intended to provide a more accurate result 
than the proxy provided by Eq. ([23]) . Values of r\ can 
be calculated for RGB stars at different evolution stages 
(Appendix [Bj}. They show that &v mt is less dominated by 
the core rotation for early RGB stars compared to more 
evolved stages. 

As the value of rj results from a balance between the 
contributions of p and g modes, we can assume that its 
value depends on the mean pressure and gravity radial or- 
ders. We therefore propose a phenomenological proxy of 
the parameter rj, justified by a simple model presented in 
Appendix IB.21 as a function of global seismic parameters. 
The analysis of the integrated kernels calculated for red 
giant models at different evolutionary stages provides the 
fit 



j/^Arii 

600 r/~l+ 7 ^ - 

where Alii is the period spacing of gravity modes de- 
rived from the fit of the mixed modes and 7 ~ 0.65 (see 
Appendix IB.2j) . The model indicates that the coefficient r\ 
decreases towards unity when a star evolves along the RGB, 

605 due to the significan t increase of the gravity radial orders 
(jMosser et alJl2012dl . The model, based on RGB stars, is 
assumed to be valid also for clump stars, since it only relies 
on global properties of the oscillation eigenfunction. 

In all cases, the maximal splitting of g-m modes is 

610 highly dominated by the core rotation, and its measure 
is close to the c ore ro tation. Furthermore, as found by 
iBeck et all (|2012h and iDeheuvels et al.l lj2012ah . we note 
that there is no immediate link between the minimum split- 
ting (1 — X)Su I0 t of p-m modes (Eq. ((!])) and the surface 

615 rotation. The minimum splitting measured for p-m modes 
is still strongly dominated by the core rotation. Extracting 
the surface rotation, which is supposed to be small, requires 
a very accurate description of the kernels, which is out of 
the scope of this work dedicated to ensemble asteroseismol- 

620 ogy 



6. Discussion 

The relation between the maximum rotational splitting and 
the mean core rotation period (Eq. ((23)) ) allows us to re- 
visit the scaling relations established in Sectional We have 
625 to reiterate that Eq. (fM| is based on a strong hypothesis, 
resulting from the linear relation between the splitting and 
the core rotation rate. A high r adial different i al rot ation 
profile in the core, a s shown by iGouoil et al.l <|2012T) and 
iMaraues et al.l (|2012r i. would invalidate the relation. 



630 6.1. Internal angular momentum transfer 

The scaling relations in Eqs. (j6]) and (JT]) can be written in 
terms of (T rot ) c rather than 6v I0 t- We find, for RGB and 
clump stars, respectively: 



(T rot ) c oc R°- 7±0 - 3 (RGB) 
(T rot ) c oc R 1A±0A (clump) 
(T rot ) c oc i? 15±0 - 4 (2nd clump). 



(26) 
(27) 
(28) 



from the fact that the correction factor 77 of Eq. (j25|) is 
close to unity and does not show important variation. As 
a consequence, the regime seen in the rotational splitting 
translates into a similar regime of the mean core rotation 
period (Fig. 640 

6.1.1. On the RGB 

On the RGB, the spinning down of the core is moder- 
ate, much smaller than a variation in R 2 expected in case 
of homologous spinning down at constant total angular 
momentum. Angular momentum is certainly transferred 645 
from the core to the envelope, in order to spin down the 
core. However, a strong differential rotation profile takes 
place when giants ascend the RGB (|Marques et al.l [20121 : 
iGoupil et al.ll2012ft . 



(25) 6.1.2. Clump stars 



650 



The extrapolation of the fit reported by Eq. (|2"6l to a typ- 
ical stellar radius at the red clump shows that cores of 
clump stars are rotating six times slower. This slower ro- 
tation can be partly explained by the core radius change 
occurring when helium fusion ignition removes the degen- 655 
eracy in the core. This change, estimated to be less than 
50% (|Sills fc Fmsonncau itl i2ooor> . can however not be re- 
sponsible for an increase of the mean core rotation period 
as large as a factor of six. As a consequence, the slower rota- 
tion observed in clump stars indicates that internal angular 660 
momentum has been transferred from the rapidly rotating 
core to the slowly rotating envelope (Fig. [HJ) - 

6.1.3. Comparison with modeling 



We note that the absolute values of the exponents are sim- 
635 ilar to the exponents found in Eqs. ([6]) and ([7]). This comes 



T he comparison with modelin g reinforces this view (Fig. 1 
of ISills fc Pinsonneaultl l2000l ). Their evolution model as- 
sumes a local conservation of angular momentum in radia- 
tive regions and solid-body rotation in convective regions. 
It provides values for the core rotation in a 0.8 M@ star of 
about 50 days on the main sequence, about 2 days on the 
RGB at the position of maximum convection zone depth in 
mass, and about 7 days in the clump. This means that, even 
in a case where the initial rotation on the main-sequence 
is slow (certainly much slower than the main-sequence pro- 
genitors of the red giants studied here) and where angular 
momentum is massively transferred in order to insure that 
convective regions rotate rigidly, the predicted core rotation 
periods are much smaller than observed. The expansion 
of the convective envelope provides favorable conditions 
for internal gravity waves to transfer internal momentum 
from the core to the envelope to sp i n down the core rota- 
tion (IZahn et al.lll997t lMathisll2009l ). iTalon fc Charbonnell 
(|2008 ) have shown that the conditions are favorable for 
these waves to operate at the end of the subgiant branch 
and during the early- AGB phase. There is observational ev- 
idence that the spinning down should have been boosted in 
the upper RGB too. 

The comparison of the core rotation evolution on the 
RGB and in the clump shows that the angular momentum 
transfer is not enough for erasing the differential rotation 
in clump stars. The line representing an evolution of (T ro t)c 
with R 2 extrapolated to typical main-sequence stellar radii 
gives a much more rapid core rotation than the extrapola- 
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Fig. 9. Mean period of core rotation as a function of the asteroseismic stellar radius, in log-log scale. Same symbols and 
color code as in Fig. [HI The dotted line indicates a rotation period varying as R 2 . The dashed (dot-dashed, triple-dot- 
dashed) line indicates the fit of RGB (clump, secondary clump) core rotation period. The rectangles in the right side 
indicate the typical error boxes, as a function of the rotation period. 



tion from the RGB fit. This indicates that the interior struc- 
ture of a red-clump star has to sustain, despite the spinning 

695 down of the core rotation, a significant differential rotation. 
This conclusion, implicitly based on the assumption of to- 
tal angular momentum conservation, is reinforced in case 
of total spinning down at the tip of the RGB. However, 
the large similarities of the values of the core rotation pe- 

700 riod observed in clump stars, together with an evolution of 
{r ro t) c close to R 2 (Eq. (|2"?| ), should imply that a regime is 
found with a core rotation of clump stars much more rapid 
than the envelope rotation but closely linked to it. 

6.2. Mass dependence 

705 Wc have calculated, for different mass ranges [Mi , Ma] , a 
mean core rotation period defined by 

E M 2 Irp V p — r 
, Mi \ J rot/c Jt 

{{J-Tot)c)[M lt M2} - „ Ma D _ r W 

where r is the exponent given by Eqs. (|26l) or (l27|) . depend- 
ing on the evolutionary status. 

710 This expression allows us to derive a mean value even 
for RGB stars, in the mass range [1.2, 1.5M ] where the 
RGB star sample can be considered as unbiased. We do not 
detect any mass dependence. The situation changes dras- 
tically for clump stars, with a clear mass dependence: the 

715 mean value of (T rot ) c is divided by a factor of about 1.7 
from 1 to 2 Mq . This reinforces the view that angular mo- 
mentum is certainly exchanged in the upper RGB, since it 



may indicate a link between the amount of specific angular 
momentum transferred and the evolution time: high-mass 
red giants evolve more rapidly than lower mass stars, loose 720 
less mass, and keep a more rapidly rotating core after the 
tip. 

7. Conclusion 

Rotational splittings were measured in about 300 red giants 
observed dur ing more than two years with Kepler. As first 725 
measured bv iBeck et alJ (|2012f l for three red giants in the 
early stages of the RGB, a strong differential rotation is 
noted for all these red giants. 

We have first shown that the rotational splitting pat- 
tern, modelled as a function of the mode frequency, is 730 
largely independent of the stellar evolutio n. The empiri- 
cal pattern found by iMosser et al.l (|2012bh has been used 
and verified in a large set of stars and has proven to be 
very efficient for analysing the splittings. Independent of 
any modeling, we have shown that the scaling relations ob- 735 
served for the maximum rotational splittings in RGB stars 
may suggest that transfer of angular momentum must oc- 
cur in their interiors. 

Then, assuming that the relation between the rotational 
splitting and the rotation rate is linear, we have shown that 740 
the measured splittings provide an estimate of a rotation 
period representative of the mean core rotation. We observe 
that this period is larger for clump stars compared to the 
RGB. This requires a transfer of angular momentum in the 
star to spin down the core. Despite the angular momentum 745 
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loss expected at the tip of the RGB, the core rotates more 
rapidly in clump stars than expected from an evolution as 
the square of the radius. This indicates a strong differen- 
tial rotation in clump stars as well as in RGB stars. In 
750 other words, the mechanism responsible for the redistribu- 
tion of angular momentum is efficient enough to spin down 
the mean core rotation but with a time scale too long for 
reaching a solid rotation. 

The indirect estimate of the specific angular momentum 
755 shows that massive red giants observed in the secondary 
clump have a significantly higher specific angular momen- 
tum than in the main red clump. 

This ensemble asteroseismic analysis of rotation in red 
giants will have to be extended to subgiants, since sub- 

760 giants also show mixed modes that give access to the inner 
rotation profile. As Kepler continues to observe, we will 
have access to longer observation runs. This will provide 
more resolved observations of the rotational splittings at 
low frequency, so we hope to measure the mean core rota- 

765 tion on the upper part of the RGB and, if mixed modes 
are also present, on the asymptotic giant branch. Our find- 
ings provide strong motivation for further stellar modeling 
including rotation. 



Appendix A: How rotational splittings are fitted 

770 A.l. Large separation and gravity mode spacing 

The first step for identifying the red giant oscillation spec- 
trum is, as for all stars showing solar-like oscillations, the 
correct identification of the radial mode pattern, in or- 
der to locate precisely the location of the theoretical pure 

775 dipole pressure modes. The fit of the radial modes depends 
mainly on the accurate determination of the large separa- 
tion. According to the un iversal red giant oscillation pat- 
tern (jMosser et al.ll20lTbl) . the surface offset and the cur- 
vature of the ridge are functions of the l arge separation. In 

780 practice, small residuals due to glitches (Miglio e t al.ll20To[ ) 
can induce a frequency offset of about, typically, Aj//50. 
Thus, a second free parameter, simply a frequency offset, 
or equivalently an offset of e less than 0.02 (Eq. (fT])), is 
useful for providing the best fit of the radial ridge. The lo- 

785 cation of the dipole ridge with respect to the radial ridge 
is given by the small separatio n 01 (Eq. (fTj)), which is a 
function of the large separation (jMosser et al.ll2011bh . 

The fit of the mixed-mode pattern is based on two 
free parameters: the period spacing Alii and the couplin g 

790 constant q 7 as defined by Eq. 9 of (jMosser et all 120123) . 
wh ich closely f o llows the formalism of mixed modes given 
by lUnno et al.l (jl989). In order to determine Alii on the 
RGB, it is worthwhile to consider that this period is a func- 
tion of the large separation. For the low-mass stars of the 

795 RGB with a degenerate helium core, a convenient proxy is 
given by the polynomial development 

Alii = 62.5 + 1.40 Av + 0.081 Av 2 (A.l) 

wi th Av in nHz and AIT in s, according to Fig. 3 
of M 

osser et ail (j2012cD . When the rotational splitting is 
800 larger than half the mixed mode spacing at f max , this step 
cannot be done independent of the next one. 
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Fig. A.l. Fit of rotational splittings, for the RGB star KIC 
6144777, with an echelle diagram as a function of the re- 
duced frequency v j Av — (n p + e). The radial orders are 
indicated on the y-axis. Radial modes (highlighted in red) 
are centered on 0, quadrupole modes (highlighted in green), 
near —0.12 (with a radial order n p — 1), and I = 3 modes, 
sometimes observed, (highlighted in hell blue) near 0.20. 
Rotational splittings are identified with the frequency of 
the to — component given by the asymptotic relation 
of mixed modes, in /iHz. The fit is based on peaks show- 
ing a height larger than eight times the mean background 
value (grey dashed lines). In order to enhance the appear- 
ance of the multiplets, highest peaks have been truncated; 
to enhance the short-lived radial and quadrupole modes, 
a smoothed spectrum is also shown, superimposed on the 
corresponding peaks. 



A. 2. Rotational splittings 

Great care must be taken to disentangle the splittings from 
the mixed mode spacings. Three major cases have to be 
considered for fitting the rotational splittings. 805 

- If splittings are small and almost uniform with fre- 
quency, except the modulation depicted by 1Z (Eq. |4]) , then 
the estimate is straightforward. The unknown stellar incli- 
nation can be derived from the mode visibility, which de- 
pends on the azimuthal order to. According to the probabil- 810 
ity of having an inclination i proportional to sini, in most 
cases doublets with to = ±1 are observed. Note that, even 

if the components to = — 1 and +1 have the same visibility, 
they may in practice present different heights, due to the 
stochastic excitation of the modes. Such splittings smaller 815 
than the mixed-mode spacings are seen in the lower stages 
of the RGB and in the clump (Fig. lA.lj) . 

- If apparent splittings at v max seem to increase with 
increasing frequency, then the most plausible solution is 
that Sv^ot is close to half the mixed-mode spacing at ^ max . 820 
These apparent splittings result in fact from a mixing of 
the splittings embedded with the spacings. Such a situation 
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Fig. A. 6. Gravity echeile diagrams of the two RGB stars KIC 5858947 and 11550492. The x-axis is the period l/v 
modulo the gravity spacing Alii; for clarity, the range has been extended from —0.5 to 1.5 Alii. The size of the selected 
observed mixed modes (red diamonds) indicates their height. Plusses give the expected location of the mixed modes, 
with 772 = — 1 in light blue, m = in green and m = +1 in dark blue. 



occurs when the apparent splittings are composed of the 
777 = ±1 component of the mixed mode order 7i m and of the 
825 777 = =pl component of the adjacent orders 77 m ± 1. The true 
splittings, significantly larger than the apparent splittings, 
are almost uniform for g-m modes. This uniformity is used 
for iterating the solution. Such cases occur most often for 
RGB stars with Av in the range [9 - 12/jHz] (Fig. lA.4|) . 



830 - If apparent splittings seem very irregular, then the 
most plausible solution is that <5^ ro t is much larger than 
half the mixed- mode spacing at i/ max . In fact, the apparent 
splittings are complex structures resulting from a mixing of 
components of two or three different mixed-mode orders. A 

835 careful visual inspection is necessary to disentangle them. 
The mixed-mode asymptotic expression and the empirical 
expression of the rotational splitting are accurate enough 
for resolving complex cases that occur for RGB stars with 
Av < 9/jHz (Fig. Eg}. 



840 We have us ed gravity eche l le dia g rams to represent th e 
mixed modes ([Bedding et all 120111: iMosser et all l2012cl) . 
Due to the complexity of the features caused by embedded 
splittings and mixed modes spacings, the echelle diagrams 
cannot be used to identify the rotational splittings, but are 

845 useful for improving the accuracy of the fit. In the examples 
shown (Fig. IA.6I) . a 10-s shift between the periods of the 
observed and modeled peaks correspond to an accuracy in 
frequency of about Ai//100. 



A. 3. A dipole mode forest? 

The complete fit of the rotational splittings is based on 850 
three parameters: the maximum splitting 5v m t and the two 
parameters A and (3 entering the definition of 1Z. The best 
fit is provided by correlating the observed multiplets with 
synthetic multiplets. 

Since the parameters A and /3 are found to vary in nar- 855 
row ranges, the solution for inferring 8v vot (and simultane- 
ously Alii on the RGB with low Av) is based on consider- 
ing them as constants. As a result, five free parameters are 
enough for fitting the whole red giant oscillation spectrum. 
Variation of A and (3 allows a better fit. The stellar inch- 860 
nation can be derived from the ratio of the visibility of the 
777 = ±1 components compared to the central component. 

In a typical spectrum, more than 30 mixed-mode orders, 
representing about 60 to 120 individual modes with a height 
larger than eight times the background are simultaneously 865 
fitted. The typical accuracy of the fit, of about Aj//200 or 
better, is enough for avoiding any confusion in almost all 
cases, except for the most evolved RGB stars. 

Finally, with the identification of the mixed mode spac- 
ings and of the rotational splittings, the dipole mode forest 870 
becomes a well-organized garden a la frangaise. 

Appendix B: Two-layer model 

B. l. Core and surface contributions 

In order to estimate the contribution of the core and surface 
rotation, we simplify the stellar stratification to a 2-layer 875 
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Fig. A.2. Same as Fig. [O] for the RGB star 5858947. In 12 
such a spectrum where the total splitting 28v m t is equal 
to half the mixed- mode spacing at v max , the fit allows to 
correctly identify the multiplets. 

4770846 - Ay = 5.48 ) LiHz - An, = 308.1s - v Tol = 65nHz 





0.2 0.4 

v / hv - (n p +E) 



0.2 0.4 
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Fig. A. 3. Same as Fig. IA.11 for the clump star KIC 
4770846. The apparent low quality of the fit for p-m 
modes at large freq uency is due to their short lifetimes 
(jBaudin et al.l 



model. We denote by 8v c and 8v s the rotational frequency 
of the core and at the surface, respectively, and 8v g /2 and 
Svp the measured splitting on g and p modes, respectively. 



Fig. A.4. Same as Fig. ED for the RGB stars KIC 9267654 
and KIC 10866415, where the total splitting 28v rot is nearly 
equal to the mixed-mode spacing at ^ max . Apparent narrow 
multiplets are artifacts due to close combinations between 
components of different mixed-modes radial orders. 



The factor 1/2 in 8v g /2 accounts for the Ledoux coefficient. 
The contributions of the surface and of the core are written 



8v„ 



x g 8v c + (1 — x g ) 8v a 



8v v = x p 8v c + (1 - x p ) 8u a 



(B.l) 



The coefficient a: p and derived from the rotational 

kernels. From the solution 



8v c 



1 Xp Xp 

—8v s + 



1 



-8v n 



(B.2) 
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11550492 - Aw = 8.70 J uHz - An i = 74 



=390nHz 




v / Av 



Fig.A.5. Same as Fig. EH for the RGB star KIC 
11550492. The non-negligible amplitudes of the m = com- 
ponents complicate the analysis. 



7 g | are constant not so far from unity. Hence, neglecting 
Sv s in Eq. (jB.ll) . we derive: 



915 



Sv c 



Tp^p 
7g n g 



1-7- 



Sv e 



(B.5) 



The radial order n p and n g have to be estimated at the 
frequency v max where the oscillation amplitude is maxi- 
mum. Then, we can derive that r/ = 5v c /6v g is related to 
the global seismic parameters Av and Alii, where Ani is 
period spacing of gravity modes: 



920 



7] ~ 1 + 7 



.All! 



Av 



(B.6) 



The fit of the integrated kernels calculated at different evo- 
lutionary stages gives 7 ~ 0.65. This phenomenological re- 
sult based on a simple two-layer model is to be considered 
as a proxy only. 
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885 and from the observation of the splitting of p-m modes 
indicating 5v p ~ Sv g /A (a factor of about 1/2 comes from 
1 — A in Eq. Q, an another factor of 1/2 comes from the 
Ledoux coefficient), one derives that the measure of 5v g is 
an indicator of the core rotation 

890 Sv c = T] Sv g . (B.3) 

For an RGB star at the bump with Av — 5 /iHz, the values 
x p and x g derived from the kernels give r\ = 1.06 ± 0.04, 
ver y close to unity. A less evolved star, as considered 
bv lBeTk et al.l (|2012t ). whose mixed modes correspond to 

895 much smaller r a dial gr avity orders, has 7/ = 1. 45tg. 15- 
iDeheuvels et al.l (|2012al) derived a similar result for a gi- 
ant with Av ~ 29 /iHz at the bottom of the RGB. This 
shows that Sv TOt is less dominated by the core rotation for 
early RGB stars. One also derives that, in all cases, the sur- 

900 face rotation Sv s is small, and that measuring it precisely 
from the g-m mode-splitting is not possible. 

B.2. Link to the eigenfunction properties 

The value of the coefficients x g and x p introduced in 
Eq. (|B.1[) can be approximated by the expression of the 

905 rotational splitting (Eqs. © and ©)• Basically, the inte- 
gration of the wave function has a contribution varying as 
the number of nodes in the core and in the envelope, respec- 
tively. As a consequence, x g oc n g and x p oc n p . In order 
to more precisely account for the complex form of the wave 

910 function, we suppose: 

f x g = 7g n g / ( 7g n g + 7 p n p ) , B 4 , 

\ x p = 7 p n p / (7 g n g + 7 p n p ) 

with 7 g < to account for the negative value of n g . The 
validity of this development implicitly assumes that 7 P and 
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